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A.  Introduction 


Water  ice  clouds  form  in  the  high  latitude  summer  mesosphere  at  altitudes 
between  80  and  87  km  when  the  temperature  falls  below  150  K  [Gadsden  and 
Schroder,  1989;  Thomas,  1991],  When  the  ice  particles  that  make  up  these  clouds 
grow  large  enough  to  scatter  light  efficiently,  which  usually  happens  at  the  base  of  the 
cloud  layer  (80  -  84  km),  the  clouds  become  visible  to  ground-based  observers  during 
twilight  when  the  lower  atmosphere  is  in  shadow.  These  spectacular  “night-shining” 
clouds  are  then  termed  noctilucent  clouds  (NLCs).  They  are  also  observed  from  space 
by  limb  scanning  satellites  when  they  are  referred  to  as  polar  mesospheric  clouds 
(PMCs). 

The  first  reported  sightings  of  NLCs  were  in  June,  1885,  when  it  was  quickly 
realised  that  the  clouds  occurred  at  the  astonishing  altitude  of  about  83  km.  There  is 
evidence  of  a  secular  increase  in  the  frequency  of  their  occurrence,  and  NLCs  have 
also  started  to  appear  at  much  lower  latitudes  (40°  N  in  1 999).  This  has  naturally  led 
to  speculation  that  they  provide  an  early  warning  of  climate  change  [Gadsden  and 
Schroder,  1989;  Thomas,  1991;  von  Zahn,  2003],  NLCs  are  composed  of  ice 
particles,  and  since  they  occur  in  an  extremely  dry  region  of  the  atmosphere,  their 
changing  climatology  must  indicate  either  that  the  mesosphere  is  cooling  or  that  the 
concentration  of  H2O  is  rising  (probably  because  of  increased  CH4  emissions). 

The  clouds  tend  to  form 


Concentration  /  cm" 


around  the  cold  mesopause 
(«  87  km),  where  the 

condensation  nuclei  are  either 
D  region  proton-hydrates  or 
metal-containing  “smoke 
particles”  produced  from 
meteoric  ablation.  As  the 
particles  grow  they  settle 
gravitationally,  eventually 
sublimating  in  the  warmer 
mesosphere  below  82  km. 
Only  the  largest  ice  particles  (> 
20  nm)  give  rise  to  the 
phenomenon  of  NLC  and 
PMC.  However,  there  is 
thought  to  be  a  much  larger 
population  of  sub-visible  small 
particles  (<  20  nm)  situated 
above  the  visible  layer  which 


Figure  1.  Annual  average  layers  of  Fe,  Na,  K  and  Ca,  PMC.  However,  there  is 
measured  at  several  lidar  observatories  in  the  U.S.  and  thought  to  be  a  much  larger 
Europe  [Plane,  2003],  The  height  range  in  which  population  of  sub-visible  small 

PMSE/NLC  nucleate,  sediment  and  sublimate  is  shown  particles  (<  20  nm)  situated 
for  comparison.  above  the  visible  layer  which 

give  rise  to  strong  VHF  radar  backscatter  signals,  termed  Polar  Mesospheric  Summer 
Echoes  (PMSEs).  PMSE  and  NLC  seem  to  be  a  closely  related  phenomenon  [von 
Zahn  and  Berger,  2003],  The  refractive  index  of  a  mesospheric  radar  is  determined 
by  sharp  gradients  in  the  free  electron  density,  which  implies  that  the  particles  giving 
rise  to  PMSE  must  be  charged.  Indeed,  the  particles  often  become  negatively  charged 
because  of  electron  capture.  However,  recent  measurements  with  rocket-bome 
particle  detectors  have  revealed  pronounced  layers  of  positively-charged  particles. 


together  with  large  enhancements  of  electrons  [Blix,  1999],  The  clear  implication  of 
these  observations  is  that  there  can  be  very  efficient  photoemission  of  electrons  from 


the  ice  particles  [Rapp  and  Lubken,  1999].  Furthermore,  since  pure  ice  has  a  high 
work  function  (~  8  eV),  photoelectric  emission  must  result  from  metal  atoms  such  as 
Na  and  K  that  are  embedded  in  the  NLC  particles. 

NLC  and  PMSE  occur  in  a  region  of  the  upper  mesosphere  where  there  are 
layers  of  metal  atoms  (Figure  1),  produced  by  the  ablation  of  the  roughly  50  -  120 
tonnes  of  interplanetary  dust  that  enters  the  atmosphere  daily  [Plane,  2003],  There  is 
growing  evidence  that  mesospheric  clouds  interact  with  the  meteoric  metal  layers, 
with  the  result  that  the  metal  layers  are  significantly  depleted  in  the  high  latitude 
summer  mesosphere.  Simultaneous  measurements  of  Fe  and  NLC  at  South  Pole  with 
a  two  channel  lidar  revealed  dramatic  bite-outs  in  the  Fe  layer  coinciding  with  NLC 
[Plane,  2004],  Measurements  of  die  Na  layer  at  Andoya,  Norway  (69°N)  (Jo  She, 
University  of  Colorado,  using  the  WEBER  lidar  funded  by  AFOSR),  at  South  Pole 
and  at  Sondrestrom,  Greenland  (67°N)  have  shown  that  the  Na  layer  is  also 
substantially  depleted  throughout  the  summer  period  below  90  km  -  the  region  in 
which  ice  particles  exist.  Finally,  measurements  of  the  K  layer  above  Spitsburgen 
(78°N)  show  that  the  metal  is  also  substantially  depleted  below  90  km,  and  the  onset 
of  low  K  densities  corresponds  to  the  onset  of  PMSE  [Hoffher  et  al.,  2003], 

One  possible  mechanism  for  NLC-metal  layer  interaction  is  the  direct  uptake 
of  the  metal  atoms  and  other  metallic  species  on  the  surface  of  the  ice  particles. 
However,  uptake  on  ice  particles  must  compete  with  the  input  of  fresh  metal  atoms 
from  meteoric  ablation,  as  well  as  vertical  transport  of  metals  into  the  ice  cloud.  The 
first-order  heterogeneous  loss  of  a  species  on  a  surface  (feet)  is  described  by 

Kt=i-r-c-A  (1) 

where  y  is  the  uptake  coefficient,  defined  as  the  probability  of  permanent  removal 
from  the  gas  phase  upon  collision  with  the  surface,  and  is  measured  under  appropriate 
conditions  in  the  laboratory;  c  is  the  root-mean-square  velocity  of  the  species  in  the 
gas  phase  and  ,4  is  the  is  the  reactive  volumetric  surface  area.  Taking  the  example  of 
atomic  Na:  the  average  meteoric  Na  input  flux  is  about  3  x  103  atoms  cm'2  s’1  [Plane, 
2004],  if  this  input  is  spread  over  40  km  (the  typical  vertical  range  of  ablation)  the 
input  rate  is  8  x  10"4  atoms  cm'3  s'1.  To  exceed  this  rate,  if  A  =  5  x  10'7  cm2  cm'3  and 
T  =  140  K  (typical  of  the  summer  mesosphere[fierger  and  von  Zahn,  2002]),  then  y 
needs  to  be  larger  than  about  0. 1 . 

The  alkali  metals  undergo  a  sequence  of  reactions  with  O3,  H2O  and  CO2  in 
the  upper  mesosphere  to  produce  the  relatively  stable  bicarbonate  species,  NaHCCb 
and  KHCO3  [Plane,  2004],  These  have  very  large  dipole  moments  (>  6  Debye),  and 
so  are  almost  certainly  lost  with  y  of  unity  on  the  polar  surface  of  ice  particles. 
However,  there  is  no  a  priori  reason  to  suppose  that  the  neutral  metal  atoms  will  be 
lost  upon  the  ice  surface  with  such  high  efficiency.  Indeed,  a  previous  study  of  the 
uptake  of  atomic  O  on  ice  films  shows  that  y  is  only  2  x  10"5  under  mesospheric 
conditions  [Murray  and  Plane,  2003a;  Murray  and  Plane,  2003b].  Hence 
measurement  of  the  uptake  coefficients  of  Na  and  K,  under  conditions  pertinent  to  the 
polar  summer  mesosphere,  are  required. 
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B.  Description  of  the  Work 

Laboratory  Measurements 


The  liquids  cooled  fast  flow  tube  used  to  measure  the  uptake  of  Na  and  K 
atoms  has  been  described  in  detail  elsewhere  [Murray  and  Plane,  2003b],  The  flow 
tube  is  constructed  from  a  borosilicate  glass  tube  of  inner  diameter  1 .98  cm,  which 
can  be  cooled  to  between  80  and  210  K  with  temperature  gradients  of  less  than  10  K. 
The  cooled  section  is  about  70  cm  long.  At  the  downstream  end  of  the  tube  there  is  a 
cross  piece  where  the  metallic  species  were  observed  by  resonance  fluorescence.  The 
alkali  atoms  were  generated  in  the  gas  phase  at  the  upstream  end  of  the  tube  and 
entrained  in  a  flow  of  bath  gas  (N2  or  He)  that  transported  them  down  the  tube.  Ice 
films  were  prepared  by  depositing  water  vapour  from  the  gas  phase  onto  the  cooled 
walls  of  the  flow  tube.  Water  vapour  was  introduced  into  the  flow  tube  via  a  heated 
injector  which  was  retracted  along  the  axis  of  the  tube  at  a  constant  rate,  in  order  to 
deposit  an  even  ice  film  thickness  along  the  length  of  the  tube.  These  ice  films  were 
deposited  at  90  K  where  they  were  most  likely  to  be  of  the  amorphous  phase,  and 
could  then  be  annealed  at  higher  temperature  in  order  to  control  the  phase  and 
morphology  of  the  ice  film.  The  ice  films  were  extensively  characterised  using  a  Kr 
gas  adsorption  isotherm  technique  and  applying  the  BET  analysis  [Brunauer  et  al, 
1938],  These  measurements  showed  that  the  surface  area  of  freshly  deposited  ice  at 
90  K  was  around  2500  times  larger  than  the  geometric  surface  area,  and  that  this 
surface  area  enhancement  decreased  to  just  90  if  the  ice  films  were  annealed  at  90  K 
[Murray  and  Plane,  2003b].  Measurements  of  the  uptake  of  the  metals  were  made  on 
ice  films  of  both  amorphous  and  cubic  crystalline  phase,  since  these  are  the  phases  of 
ice  thought  to  be  most  important  in  the  summer  mesosphere. 


Figure  2.  A  schematic  drawing  (approximately  to  scale)  of  the  fast  flow  tube  employed  to  study 
the  uptake  of  atomic  K  and  Na  on  low  temperature  ice.  (PMT  =  photomultiplier). 
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Gas-phase  Na  (or  K)  atoms  were  generated  using  an  oven  which  was  heated  to 
around  620  K  (or  520  K),  where  the  vapour  pressure  of  the  alkali  metals  is  significant 
(«  1  x  1015  atom  cm'3).  The  oven  consisted  of  a  stainless  steel  tube  (19  mm  diameter), 
wrapped  with  heating  tape.  The  oven  was  then  welded  to  a  6  mm  diameter  tube  of 
500  mm  length,  which  served  as  an  injector  (Figure  2).  This  injector  was  heated  to 
670  K  (or  570  K)  using  a  heating  element  wound  around  the  outside  of  the  6  mm 
tube;  the  heater  was  in  turn  encased  in  a  double-layered  steel  jacket  through  which 
chilled  cooling  water  flowed.  This  design  ensured  that  the  central  6mm  tube  through 
which  the  metal  atoms  flowed  was  sufficiently  hot  for  the  atoms  to  remain  in  the 
vapour  phase,  while  the  outside  skin  of  the  injector  was  maintained  slightly  below 
room  temperature. 

The  alkali  metals  were  therefore  injected  directly  into  the  cooled  section  of  the 
flow  tube,  where  they  were  entrained  in  the  large  flow  of  carrier  gas  (N2  or  He)  about 
45  cm  upstream  of  the  cross  piece.  The  relative  alkali  metal  concentration  at  the  cross 
piece  was  measured  by  resonance  fluorescence,  using  a  hollow-cathode  lamp 
operating  at  589  nm  (or  760  nm),  a  narrow  bandpass  interference  filter  and  photon 
counting. 


Figure  3.  The  first  order  loss  of  atomic  Na  (triangle),  K  (squares)  and  Fe  (circles),  in  He  (open 
symbols)  and  N2  (filled  symbols)  carrier  gases  on  the  uncoated  Pyrex  flow-tube  walls.  The  Na  and  K 
concentration  were  proportional  to  the  resonance  fluorescence  signal  measured  by  photon  counting. 
The  plots  have  been  offset  on  the  arbitrary  concentration  scale  for  clarity.  The  pressure  was  1.5  Torr 
and  the  carrier  gas  velocities  ranged  between  7  and  60  m  s'1. 


Results  and  Discussion 

The  loss  of  Na  or  K  through  uptake  on  the  flow  tube  walls  was  determined  by 
measuring  the  relative  concentration  of  the  metal  at  the  downstream  cross  piece,  while 
varying  the  flow  time  0am  (the  time  the  bath  gas  takes  to  between  the  injector  and  the 
cross  piece)  at  constant  pressure  in  the  flow  tube.  t\am  was  calculated  from  the  known 
mass  flow  rates,  pressure  and  flow  tube  dimensions,  and  corrected  for  the  interaction 
between  the  flow  dynamics  and  reactive  uptake  at  the  walls.  This  correction  arises 
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when  the  removal  of  a  gas-phase  species  on  the  walls  is  efficient,  so  that  the  loss  to 
the  walls  becomes  limited  by  diffusion  and  a  radial  gradient  forms  across  the  tube 
with  the  maximum  concentration  in  the  centre  of  the  tube.  The  centre  is  also  where 
the  velocity  of  the  bath  gas  is  greatest  due  to  the  laminar  flow  regime  (which  existed 
under  all  flow  regimes  used  here).  The  combination  of  laminar  flow  and  radial 
concentration  gradients  results  in  the  metal  species  travelling  down  the  flow  tube  at  a 
greater  velocity  than  the  plug  (or  average)  flow  velocity,  v,  of  the  bath  gas:  when 
there  is  rapid  wall  loss  then  Mam  =  1.60  k,  where  Mam  is  the  average  velocity  of  the 
metal  atoms. 

As  shown  in  Figure  3,  the  loss  of  Na  and  K  to  the  flow  tube  walls  occurs  with 
first-order  kinetics  in  both  He  and  N2  bath  gas.  The  slopes  of  the  lines  in  this  Figure 
yield  the  first-order  rate  coefficients,  &n3  and  km.  Note  that  sets  of  data  for  atomic  Fe 
removal  are  included  in  Figure  3  for  comparison.  The  first-order  loss  coefficients  for 
uptake  of  atomic  Na,  K  and  Fe  on  the  room  temperature  Pyrex  walls  of  the  flow  tube 
are  inversely  proportional  to  the  total  pressure  P,  as  shown  in  Figure  4.  The  absence 
of  curvature  in  these  plots  shows  that  the  uptake  of  the  metals  on  room  temperature 
Pyrex  occurs  at  close  to  unity  efficiency  (y  >  0.1).  To  illustrate  this  point,  has 
been  calculated  using  the  Brown  formalism  [fir  ow«.  1978]  for  the  hypothetical  case  of 
/Na  =  0.03,  and  is  plotted  as  a  function  of  P~l  with  a  dashed  line  in  Figure  4. 


Figure  4.  The  pressure  dependence  from  0.8  to  6  torr  of  fcMe  for  Na  (filled  diamonds),  K  (open 
triangles)  and  Fe  (open  circles),  in  He  bath  gas  at  293  K.  The  solid  lines  are  best  fits  through 
the  data.  The  error  bars  represent  the  standard  deviation  from  the  kUe  measurement  at  each 
pressure.  The  dashed  line  represents  a  hypothetical  line  where  ks„  was  calculated  using  the 
Brown  formalism  for  an  uptake  coefficient  of  Na  on  Pyrex  of  0.03. 


In  fact,  when  the  uptake  of  a  species  on  the  surface  occurs  with  unity 
efficiency  the  loss  to  the  walls  is  fully  controlled  by  diffusion.  Hence,  the  measured 
loss  rate  provides  a  measure  of  the  diffusion  coefficient,  DMe: 


5-81A. 

r2P 


(2) 
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where  r  is  the  flow  tube  radius  [Helmer  and  Plane,  1993].  In  a  laboratory  study  of 
the  uptake  of  a  species  on  a  surface  it  is  essential  to  have  an  accurate  value  for  the 
diffusion  coefficient,  especially  when  y  approaches  unity.  Hence,  in  this  study  the 
diffusion  coefficients  of  Na  and  K  were  determined  as  a  function  of  temperature  by 
measuring  the  uptake  of  the  metals  on  the  uncoated  Pyrex  surface  of  the  flow  tube, 
where  uptake  proceeds  very  close  to  unity  efficiency.  Figures  5  and  6  illustrate  the 
temperature-dependence  of  hue  for  Na  and  K  in  He,  between  80  and  295  K.  The 
temperature-dependence  is  well-fitted  by  the  form  P .  Table  1  summarises  the 
diffusion  coefficients  of  Na,  K  and  Fe  in  He  and  N2.  As  expected,  the  diffusion  of 
the  metal  atoms  through  He  is  considerably  faster,  and  so  this  bath  gas  was  employed 
for  the  uptake  measurements  on  ice.  Note  that  the  /?  values  in  Table  1  are  close  to 
the  range  of  1.7  -  1.83  expected  from  theory  [Maitland  etal.,  1981], 

Table  1.  Measured  diffusion  coefficients  of  metal  atoms  in  He  and  N2 


M 

Bath 

gas 

- K - - 

^Me  expt 

Torr  cm2  s'1 

P 

(80  -  295  K) 

Fe 

He 

366  ±  17 

1.85  ±  0.07 

Fe 

n2 

112  ±  4 

- 

Na 

He 

286  ±  13 

1.68  ±0.04 

Na 

n2 

125  ±1 

- 

K 

He 

247  + 15 

1.6910.07 

K 

n2 

88  ±2 

- 

The  uptake  rates  of  Na  and  K  on  ice  are  also  shown  as  a  function  of 
temperature  in  Figures  5  and  6.  The  ice  was  always  deposited  at  90  K  as  an 
amorphous  film.  Two  forms  of  ice  were  then  studied.  In  the  first  case,  the  film  was 
annealed  at  160  K  to  convert  it  to  cubic  crystalline  ice,  before  studying  the  uptake  of 
the  metal  atoms  at  lower  temperatures.  In  the  second,  the  ice  was  heated  from  90  K  to 
the  temperature  at  which  the  uptake  was  measured.  It  would  thus  have  been 
amorphous  up  to  the  transition  temperature  around  130  K.  The  figures  show  that 
there  is  essentially  no  difference  in  the  uptake  rate  on  the  two  types  of  ice  film.  This 
implies  that  uptake  must  be  very  efficient  -  otherwise,  the  much  greater  surface  area 
of  the  amorphous  ice  would  lead  to  significantly  larger  kMe  values  [Murray  and  Plane, 
2003b]. 

Figures  5  and  6  also  show  that  there  is  no  significant  difference  in  the  uptake 
rate  on  ice,  compared  with  the  Pyrex  surface.  From  this  observation,  and  application 
of  the  Brown  formalism  [Brown,  1978],  we  deduce  that  y(Na  on  ice)  >  0.09  and  y(K 
on  ice)  >  0.05  over  the  temperature  range  80  to  160  K.  These  lower  limits  arise 
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because  as  uptake  becomes  very  efficient,  diffusion  to  the  walls  of  tube  starts  to 
dominate  the  removal  rate  of  the  atoms  as  they  flow  down  the  tube. 

This  effect  is  illustrated  in  Figure  7,  which  shows  that  as  the  removal  rate 
constant  (fcNa)  approaches  the  diffusion-controlled  limit  (faiff),  then  yNa  approaches 
unity  very  rapidly.  For  instance,  when  k^a  is  96%  of  km,  yNa  is  still  only  about  0.1. 
Because  £Na  /  km  cannot  be  measured  with  much  greater  accuracy  than  about  4%,  yNa 
values  greater  than  about  0.1  cannot  be  measured,  only  a  lower  limit.  This  limitation 
applies  to  all  techniques  where  thermally  equilibrated  uptake  from  a  bulk  gas  is  being 
studied  (in  contrast  with  beam-scattering  techniques;  however,  these  suffer  from  the 
disadvantage  of  higher  energy  impacts  at  the  surface).  Nevertheless,  because  uptake 
for  Na  and  K  was  found  to  be  diffusion-controlled  over  a  large  temperature  range  (80 
-  160  K),  it  is  likely  that  y  is  actually  very  close  to  unity.  For  the  purpose  of 
geophysical  modelling  described  below,  we  therefore  treat  y  as  1 . 


Figure  5.  The  temperature  dependence  of  the  uptake  rate  constant  for  Na  at  fixed  temperature  (1.5 
torr)  on  Pyrex  (filled  circles),  on  ice  deposited  at  90  K  and  annealed  to  160  K  (open  circles),  and  on 
ice  annealed  at  a  temperature  no  higher  than  the  temperature  at  which  the  uptake  measurement  was 
performed.  The  dashed  line  is  a  best  fit  to  the  Pyrex  data  and  shows  that  fcNa  varies  as  7'16S1  0  04. 
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Figure  6.  The  temperature  dependence  of  the  uptake  rate  constant  for  K  on  Pyrex  and  ice.  The 
key  is  the  same  as  in  Figure  5.  kK  varies  as  TlJ±  01 . 


^Na^diff 

Figure  7.  Plot  showing  the  uptake  coefficient,  Ynb,  for  Na  atoms  on  the  flow-tube  walls, 
as  a  function  of  the  first-order  uptake  rate  k^a  normalised  to  the  diffusion-controlled 
limit  The  calculation  is  performed  using  the  Brown  formalism  [Brown,  1978] 
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Atmospheric  Modelling 


The  implications  of  these  measured  uptake  coefficients  were  then  explored 
using  the  UEA  1 -dimensional  model  of  the  mesosphere/lower  thermosphere  [Plane, 
2004].  This  model  extends  from  65  to  110  km  with  a  resolution  of  0.5  km,  and  uses 
time-implicit  integration  with  a  10  min  time  step.  The  model  contains  a  full  treatment 
of  the  odd  oxygen  and  hydrogen  chemistry  in  the  mesosphere,  and  the  neutral  and 
ion-molecule  chemistry  of  sodium  The  meteoric  ablation  profile  is  calculated  using 
the  standard  ablation  equations,  with  a  meteoroid  size  and  velocity  distribution  from 
the  Long  Duration  Exposure  Facility  experiment.  The  temperature  profiles  are  taken 
from  a  recent  South  Pole  climatology  [Pan  and  Gardner,  2003],  The  vertical  eddy 
diffusion  coefficient,  meridional  wind  and  vertical  wind  are  from  a  global  circulation 
model  (http://www.acd.ucar.edu/models/SOCRATES/).  The  size  and  number 
distribution  of  PMC  ice  particles  is  taken  from  a  recent  cloud  microphysics  model 
[von  Zahn  and  Berger,  2003], 

Figure  8  shows  the  effect  of  PMC  on  the  Na  and  Fe  layers  in  January  at  South 
Pole.  In  the  height  range  of  the  PMC  shown  by  the  red  line  (84  -  88  km),  nearly  all  of 
the  metal  atoms  are  removed  when  the  cloud  is  present.  This  is  because  the  rate  of 
uptake  on  the  ice  particles,  with  a  time  constant  of  between  1  and  2  hours,  is  much 
faster  than  the  rate  at  which  the  metals  are  replenished  in  the  gas  phase  by  meteoric 
ablation  and  vertical  transport. 

Figure  9  compares  the  model  predictions  with  a  comprehensive  set  of  lidar 
measurements  made  at  South  Pole  between  1995  and  1997  by  C.  S.  Gardner  and 
coworkers  at  the  University  of  Illinois.  The  panels  compare  the  measured  and 
modelled  seasonal  variations  in  the  layer  column  abundance,  centroid  height  and  root- 
mean-square  width.  Without  including  uptake  of  Na  on  PMCs,  the  very  dramatic 
seasonal  variations  in  the  layer  cannot  be  modelled  satisfactorily.  Figure  10  shows 
contour  plots  of  file  seasonal  variation  of  Na  (the  lidar  data  is  quite  heavily  smoothed 
because  of  the  lack  of  data  at  certain  periods  of  the  year).  The  abmpt  disappearance 
of  Na  when  PMC  are  present  during  mid-summer  (Nov-Feb)  is  clear  in  both  the 
measurements  and  the  model  predictions. 

Figure  1 1  illustrates  the  correlation  coefficients  between  the  sodium  and  iron 
densities  and  temperature,  as  a  function  of  altitude.  For  Na  in  particular,  excellent 
agreement  between  the  model  and  the  measurements  is  achieved.  Note  that  the 
correlation  coefficient  is  nearly  +1  below  95  km,  where  there  is  a  strong  correlation 
between  Na  density  and  temperature:  the  removal  of  Na  increases  substantially  when 
T  falls  below  150  K  and  PMCs  form.  In  contrast,  above  95  km  the  correlation 
coefficient  is  nearly  -1 ,  because  there  is  a  strong  anti-correlation  between  Na  and  T  as 
a  result  of  ion-molecule  chemistry:  Na+  ions  are  neutralised  to  Na  by  forming  clusters 
at  lower  temperatures  and  then  undergoing  dissociative  recombination  with  electrons. 

Finally,  this  modelling  exercise  has  shown  that  the  Na  layer  can  be  used  as  a 
sensitive  marker  of  atmospheric  dynamics.  One  important  finding  is  that  the 
downward  vertical  wind  velocity  into  the  polar  vortex  during  winter  is  only  about 
25%  of  that  predicted  by  global  circulation  models,  which  also  explains  why  these 
models  overpredict  the  temperature  in  the  wintertime  polar  mesosphere  [Pan  et  al., 
2002]. 


10 


PMC  surface  area  /  cm'2  cm'3 


0  le-8  2e-8  3e-8  4e-8  5e-8  6e-8 


Na  concentration  /  cm'3 


- Na,  no  PMC 

- Na  with  PMC 

Fe,  no  PMC 
Fe  with  PMC 
- PMC  surface  area 


Figure  8.  Modelled  vertical  profiles  of  Na  (black  lines)  and  Fe  (grey  lines),  when  PMC 
absent  (solid  lines)  or  present  (broken  lines).  The  red  line  shows  the  PMC  surface  area, 
calculated  from  a  cloud  microphysics  model  [von  Zahn  and  Berger ,  2003]. 
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Figure  9.  Comparison  of  the  column  abundance,  centroid  height  and  root-mean-square  width 
of  the  Na  layer  at  South  Pole,  measured  by  the  University  of  Illinois  lidar  and  calculated  by 
the  UEA  1-D  mesosphere-thermosphere  model. 
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Figure  10.  Contour  plots  of  the  seasonal  distribution  of  atomic  Na  at  South  Pole,  measured 
by  the  University  of  Illinois  lidar  (top  panel)  and  predicted  by  the  UEA  model  (bottom  panel). 
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Figure  11 .  Comparison  of  measured  and  modelled  vertical  profiles  of  the  correlation 
coefficient  between  the  metal  atom  density  and  temperature  at  South  Pole:  (a)  Na,  (b)  Fe. 
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